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Abstract 

We present an analytic formula for the 0(a s ) corrections to the 
decay t — ► H + b for nonzero mn and m^. The Equivalence Theorem 
is used to relate these corrections to those for the process t — > W + b. 
Contrary to previously published results we find that the effect of the 
b quark mass is negligible for a wide range of mt, mn and tan (3. 
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1 Introduction 



Many extensions of the Standard Model contain more than one Higgs doublet. 
The electroweak gauge bosons only absorb one of the charged Higgs fields, 
leaving the others as physical charged scalars, so that in a two-doublet model 
the top could decay to H + b (if m t > m H + + nit,). For certain choices of 
parameters, this process dominates over the expected t — > W + b decay, so it 
is of interest to correctly calculate its QCD corrections. 

It is well known that the QCD corrections to the decay rate of a heavy 
quark into a W are of order 10% in the Standard Model |], ||. Several groups 
have undertaken to calculate these corrections for the decay into a charged 
Higgs. The effect of the soft gluons has been calculated in ref. ||, and the 
decay t — > H + bg, with real gluons only, has been studied in ref. ||]. The 
full one-loop QCD corrections have been computed, at first by neglecting the 
mass of the 6-quark in ref. H, and then with a nonzero mb in ref. || (in the 
framework of what is called the Model I - see discussion below). The claim of 
the latter paper is that the QCD correction for m t = 150 GeV, m?, = 4.5 GeV 
and a s = 0.1 is as large as —15% in the (unphysical) limit of the massless 
Higgs. However, as will be argued in the following section, according to the 
Equivalence Theorem (see |7], |8|, |S| as well as [IU| and references therein) the 



correction in this limit should be the same as the correction to the decay 
t — > W + b, which for the above values of parameters is —8.6% |], [11]. The 
purpose of the present paper is to reevaluate the first order QCD corrections 
with non-zero 6-quark mass. We use dimensional regularization for both the 
ultraviolet and the infrared divergences, which leads to simpler algebra than 
if one assigns a finite mass to the gluon, as was done in 0, |(| [10|] . Our result 
is consistent with the Equivalence Theorem. 

In a model with two Higgs doublets and generic couplings to all the 
quarks, it is difficult to avoid flavour-changing neutral currents. We therefore 
limit ourselves to models that naturally side-step these problems by restrict- 
ing the Higgs couplings |Tj|. The first possibility is to have the doublet if 2 
coupling to all the quarks, and the H\ doublet interacting with none of them. 
The vacuum expectation value of Hi = v\ will nonetheless contribute to the 
W mass, leading to an H~tb vertex of the form 

tb H~b {m t cot (3 R — cot (3 L} t (model I) (1) 
V2m w 



2 



where H~ is the physical charged Higgs, Vtb is the '33' element of the CKM 
matrix, L and R are the chiral projection operators, and tan/3 = vi/v\ is the 
ratio of the vacuum expectation values of the two Higgs. The second possibil- 
ity is to have H 2 couple to the right-handed up-type quarks (ur, cr, t R ), and 
the Hi couple to the right-handed down-type quarks. This is what happens 
in the Minimal Supersymmetric Standard Model. It is easy to show that the 
interaction Lagrangian 



leads to the vertex 

gVtb 



#2<44Vl + Hi^ti + h-c (2) 



H b {m t cot (3 R — mi, tan/? L} t (model II) (3) 



where we have numbered the models in accordance with [J^ . 

In order to simplify the following formulas we introduce dimensionless 
parameters for the scaled masses: 

m b m H m w 

e = — , X = , w= , (4) 

m t m t m t 

and write the vertex t — > H + b as 
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where from (|l|) and (^|) 



V tb b (a + b 75 ) t. (5) 



Model I : { I = CO \ ^ ~ 6 \ Model II : { » = + 6 ^ ? 

I b = cotp(l + e) I b = cotp — etanp 

The next section contains our result, which is examined in the Discussion. 
The Appendix contains some details of the calculation. 



2 QCD Corrections 

The notation we use is similar to that used in the analysis of semileptonic 
decays 0, 0. In terms of the dimensionless parameters ([§), we define the 
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following kinematic variables: 



h = 

Ps = 

p± = 

% = 



^(i-x 2 + e 2 ) 

1 



l + x 4 + e 4_ 2 ( X 2 + e 2 + x 2 e2) 



2 

2 P_ 



W 

w± 

Y 

J- id 



W ±P 3 
in 



2 W_ 
The tree level decay rate is 

r°(f - = ^|Vk| 2 [P (a 2 + b 2 ) + e(a 2 - b 2 )]P 3 



and the 0(a s ) correction is 



rW = ^ Gf ^' 2 [(a 2 + b 2 ) G + + (a 2 - b 2 ) eG_ + abG 



with 



P H + P P 3 



f (2 - X 2 - 4X 4 + 3x 6 - 26 2 - 2e 4 + 2e 6 - 4 X 2 e 2 - 5 X V 



2 Ux , 



4x 



W + h [6 - 4 In f^) j + (l - x 2 - 2e 2 + e 4 - 3 X 2 e 2 ) , 



G = -6P P 3 lne, 
and 



Li 2 (P + ) - Li 2 (P_) - 2Li 2 ^1 - £^ 



+F p ln 



'4Pfx 



2 . - Y w In e 



+ 2F w (l- e 2 )+Ap 3 i ne (i + x2 _ e2 ) 



In the limit of the zero mass of the b quark the QCD correction becomes 



limr (1) 



a s G F ml \ V t 



67T 2 



^p^cot 2 /?(2G+ + G 



(10) 



where 
G + = 



G 



+ In x In ( 1 - X 

2 



Li 2 1 -X 



+ 



2x 2 



i-x 2 

5 



1-- X Z In 1-x 



vr 2 9 3, 
1 h - lne 

3 8 4 



- 3 -n 

2 



X 



lne, 



[ID 



and we see that the mass singularities ~ In e cancel in the expression for the 
total rate (|10|) . Our result in this limit is identical to the one obtained by 
Liu and Yao [[U]] and is in agreement with the corrected version of ref. |J. 
If we further take the limit m# — > the rate becomes: 



lim r (1) 



a s G F ml \ V t 



tb\ 



6vr 2 



V2 



cot 2 /3 - - 



7T 



(12) 



which is in agreement with the conclusion of the ref. [|l(|] as well as with our 
previous result [[13|]. 

Now we would like to compare the corrections to the decay width 
r (t — > H + b) with those to T (t — > W + b). For simplicity we now take mj = 
and cot (3=1, and examine the ratio of the first order correction to the Born 
rate: 



f H (x) 

fw(w) 



r« (t -> H+b) 



r(°) (t 
r« (t 



H+b)' 
W + b) 



r(o) (t _> w+6) 



(13) 



It has been noted in ||T0[ that in the limit of the infinite top mass these ratios 
are equal: /#(0) = /W(0). On the other hand, when mu approaches m t , we 
have: 



/// X — ► ^~ 

671 



-61n(l-x' 



V + 13 



(14) 
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By comparison with [ED] we see that: 



hm |44 = 1- (15) 
*-* 1 fw(x) 

Finally, we examine the corrections in the limiting case where the mass 
of the charged Higgs is zero but the b quark mass is finite. This is of course 
unphysical, but serves as a useful check on our equations. If we choose the 
parameters a and b from @ to correspond to the couplings of the single 
Standard Model Higgs, then the Equivalence Theorem implies that the cor- 
rections are the same as in the process t — > bW in the limit of massless W 
boson and nonzero to&. The latter can be obtained by taking the limit of the 
relevant formula |], |Tl~| : 

hm rO)(t^bw) - as G ^ m 'l^l 2 



m w ->0 247T 2 y/2 

\u 2 (e 2 )-4 + ln(e)ln(l-e 2 ) 



W) 2 (l + e 2 ) 
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- 4e 2 (7 - 5e 2 + 4e 4 ) In (e) - 8 (l - e 2 ) 3 In (l - t 2 ' 

-(l-e 2 )(-5 + 22e 2 -5e 4 )| (16) 

The couplings of the Goldstone boson charged Higgs of the Standard Model 
(longitudinal W) to t and b can easily be calculated to be those of Model I, 
with cot/3 = 1. In this case, the corrections to the decay t — ► bH + are, in 
the limit m# — > 0: 



lim ft -> bH~ 



a s G F ml \V t 



tb\ 



6tt 2 

2(1 + e 2 )G° + - Ae 2 G°_ + (1 - e 2 )G° 



(17) 



where G° are limits of corresponding functions Gi for tuh = 0: 
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Cf: 



'3 



Li 2 (e 2 ) - y + In (e) In (l - e 5 
+ g + e 2 -^)ln( e) -(l-^)ln(l- e 2 ) + ^- 4 ) 



2(1 + e 2 



Li 2 (e 



7T~ 



+ In (e) In ( 1 - e 



+ (3 - e 2 ) In (e) - 2 (l - e 2 ) In (l - e 2 ) + 2(1 - e 2 ) 



--(l-e 4 )ln(e) 



(18) 



Inserting these expressions into equation (|I7|) we obtain the same formula as 
(1161). 



3 Discussion 

In Figure 1 the ratio of the first order QCD correction to the Born rate 
for the decay t — > H + b is plotted as the function of the ratio of masses 
X = mu/mt. We have chosen the set of parameters cot j3 = 1, mt = 150 GeV 
for an easy comparison with the analogous diagram in ref. 0. It can be seen 
that the effect of the mass of the 6-quark is negligible, except in the case of 
m t — rriH ~ vn^. 

There are also logarithmic corrections (~ e 2 lne) to the decay rate in 
model II, as can be seen from Figure 2. Here we compare the branching 
ratios of the decays t — > H + b and t — > W + b, taking mt = 100 GeV so that 
this plot can be easily compared with a similar one published in ref. 
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Our graph is different from theirs in that we now include QCD corrections 
to both decay rates. These corrections modify the diagram significantly only 
in the case of large values of tan /3 in Model II. It must be noted however, 
that although the corrections are relatively large, the top decays principally 
to W + b in this region of tan f3. In model I, both a and b are proportional 
to cot (3, so the log of the branching ratio as a function of tan (3 decreases 
with a slope of -2. As can be seen from figure 2, the In e corrections cancel 
to order e 2 among G + , G_ and G . However in model II, the decay rate is a 
polynomial in tan/5 with exponents -2, and 2, and the lne does not cancel 
in the "0" and "2" terms. 
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Summary 



We have calculated an analytic expression for the 0(a s ) corrections to the 
decay t — > H + b for non-zero mn and m&. The m& = limit of our result 
is the same as in |I| and the corrected version of ||; however, our full 



expression disagrees with the mi, 7^ result presented in ||. To check the 
m b dependence of our results we have compared it with the corresponding 
formula for the decay t — > W + b. 
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A Details of the calculation 



Throughout this calculation we have used dimensional regularization, work- 
ing in d = 4 — 2u dimensions. The counterterm for the vertex function has 
been calculated according to ref. 



5A 



r 



a + b Srrit a — b 5m b 



m t 



m b 



+ 



\{z t -i) + l -{z b -i) 



a + b 5m t a — b 5m b 



m t 
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and the renormalization constants are (we take the renormalization scale 
equal to the mass of the decaying quark) : 



Z t -1 

z b -i 



m t 

5m b 

m b 



3tt 
3tt 



+ 37-31n4vr- 

3 4vr 
h 37 - 3 In — 

to e z 



(19) 
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The contribution of the virtual corrections to the total decay rate is: 



virt 



a s G F ml \ V t 



tb\ 



6tt 2 



where 



V- 



V2 



PoV+ Yp 



2x 2 



(a 2 + b 2 ) V + + (a 2 - b 2 ) eV_ + abG 



;i- e 2 ) 2 (l + e 2 )-X 4 (3 + 3e 2 -2 X 2 



,-2\2 



v + y B 



x" 



with 



V = 2P 3 
+ 2P 



-- + 2 7 + 21n^ 



Li 2 P + -Li 2 P_ 



3 + — (l-e 2 + X 2 )lne 



F 2 -21neK. 



(20) 



(21) 



- 2 M^ + 7 + ln ^ 

The real gluon radiation calculation requires an integration over three 
body phase space $(£; H, b, g). This leads to divergences due to the emission 
of soft and collinear gluons, for which we use dimensional regularization (see 



ref. |H| and references quoted therein). The phase space integration we are 
considering now is analogous to the decay t — > Wbg, for which the relevant 
integrals have been listed in [jTT]| . Adding the result to the virtual gluon 
contribution ( |20|) yields the final formula @. 

We would like to add that this calculation was greatly facilitated by the 
use of algebraic manipulation programs FORM Jl6[ and Mathematica [17] . 
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Figure 1: Ratio of the first order correction to the Born rate for tan/? = 1 
and m t = 150 GeV: m b = 4.5 GeV (Model I - solid, Model II - dotted), 
m b = (dash-dotted line). 



BR{t^bH+) 
BR(t^bW+) 



tan/? 



Figure 2: Ratio of the branching ratios of t — > M/ + and t — > bW + for 
m b = 4.5 GeV: in the Model I with and without QCD corrections (solid and 
dashed lines, undistinguishable) and in the Model II (dash-dotted and dotted 
lines, respectively). 
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